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, We investigate the polarized gluon distribution Ag(x) by a global analysis of cur- 

» I ■ rent DIS data and the 7r° data from RHIC-Spin experiments. The n° data provide 

' a strong constraint on Ag(x), so that its uncertainty is reduced. Furthermore, 

new DIS data of COMPASS and HERMES play an important role in determining 
Ag(z) at large x. 
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1. Introduction 

We have investigated the polarized parton distribution functions (polarized 
PDFs) by a global analysis with deep inelastic scattering (DIS) data. 1,2 
The polarized valence- up and -down distributions, Au v (x) and Ad v (x), axe 
determined well; however, the antiquark and gluon distributions, Aq(x) and 
Ag(x), have rather large uncertainties. It is therefore expected that Ag(x) 
can be extracted from collider data in which gluon contributions dominate. 
Fortunately, double spin asymmetry for the 7r° production is measured 
by the PHENIX collaboration at RHIC. 3 Because we are interested in a 
possible constraint on Ag(x) from the asymmetry data, an attempt is made 
to determine Ag{x) by an analysis including the new data. 
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2. Global analyses of the polarized PDFs and Results 

In this analysis, we chose the following functional form as a polarized PDF 
at the initial Q 2 (= Q 2 ,): 

A/Or, Ql) = [8x v - k{x v - x»)]f{x, Q 2 ) , (1) 

where 5, k, v, and /i are free parameters, and f(x) is the unpolarizcd PDF. 
The positivity condition |A/(x, Qo))l — f{ x iQo) 1S imposed as a constraint 
especially on the large-x behavior of the polarized PDFs. Moreover, the 
antiquark distributions are assumed to be flavor SU(3) symmetric due to 
the lack of accuracy in present experimental data for flavor separation of 
these distributions. We prepare four type distributions: Au v (x), Ad v (x), 
Aq(x) and Ag(x). Q 2 dependence of the PDFs is taken into account by 
solving the DGLAP equations numerically. The polarized PDFs, strictly 
speaking the free parameters in Eq.(l), are determined by a \ 2 analysis in 
the next-to- leading order (NLO). Uncertainties of these distributions are 
estimated by the Hessian method. More details are found in Rcf. 2 on the 
analysis method and the uncertainty estimation. 

We performed two analyses. The one is by using only the DIS 
data, which include recent ones from COMPASS, HERMES, and JLab 
experiments. 4,5 ' 6 The other is by using the 7r° asymmetry data of RUN05 
at RHIC 3 in addition to the DIS data. The total number of the experimen- 
tal data is 421, in which 413 and 8 are for the DIS and tt° data, respectively. 
The value of total x 2 (/cZ.o./.) is 358 (0.89) for the analysis with only the 
DIS data, and it is 370 (0.90) for the one with the tt° and DIS data. 

From these analyses, we obtain the first moments of Ag(x): 

AG(DIS + tt°) =0.31 ±0.32, 
AG(DIS) = 0.47 ± 1.08. 

The uncertainty is significantly reduced if the n° data are used although the 
center value is slightly varied. Next, both polarized gluon distributions are 
compared in Fig. 1. The Ag(x) changes in the region 0.03 < x < 0.5, and 
its uncertainty is much reduced because of the 7r° data. It indicates that 
the extraction of Ag(x) from the DIS data is difficult because the gluon 
contribution to g\{x) is rather small. It contributes indirectly through the 
Q 2 evolution and a higher-order correction. However, the ir° production 
data are useful to improve the precision of Ag(x) determination. 

In the analysis with the 7r° asymmetry data, there is, however, a prob- 
lem of Ag(x) sign. The polarized cross section is roughly proportional to 
square of Ag(x) because the gg scattering process dominates in the low-p^ 
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Figure 1. Comparison of the polarized gluon distributions and their uncertainties. The 
solid and dashed curves show iAj(i) from the DIS and 7r° asymmetry data and from 
only the DIS data, respectively. The shaded areas are their uncertainties. 

region. Therefore, two types of solutions are allowed: positive and negative 
Ag(x). In practice, we perform an analysis of a negative Ag(x) input as 
an initial condition. The value of the minimum \ 2 for the ir° data becomes 
11.05 which is almost the same as the one of the the positive Ag(x) input. 
Therefore, the sign cannot be determined by the \ 2 values. 

In addition, we found an interesting fact that Ag(x) becomes positive 
at large x in both cases. It is caused by the DIS data of the HERMES 
and COMPASS experiments. Figure 2 shows comparison of the AAC fit- 
ting results with the asymmetry data of deuteron target. Solid curves are 
full NLO calculation, and dotted curves are obtained by eliminating the 
NLO gluon term from gf{x). The differences between these curves indicate 
the gluon contribution as a higher-order correction. For the COMPASS 
kinematics, the differences are small because the Q 2 values are larger than 
those of the HERMES and the NLO correction is smaller. However, the Q 2 
values for the HERMES data are a few GeV 2 , so that the NLO contribution 
is rather large. The differences are significant in the region 0.02 < x < 0.1 
where deviations between fit results and data exist. In order to explain the 
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Figure 2. Comparisons of the AAC fit with the asymmetry data Af(x, Q 2 ) of the HER- 
MES and COMPASS. 2 
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Figure 3. z dependence of the polarized gluon coefficient function AC 9 (z) and gluon 
distribution Ag(x/z) in the convolution integral for g\(x). 2 

data, the NLO gluon term should be positive. 

For obtaining a positive gluon term in such an x region, Ag(x) must be 
positive at large x. The gluon term is given by the convolution integral with 
the coefficient function AC g (z): dz / zAC g (z)Ag{x / z). The behavior 

of these functions is shown in Fig. 3. The coefficient function is positive in 
the region 0.02 < z < 0.7. To obtain the positive gluon term at x = 0.05 
where the deviation from the HERES data exists, the gluon distribution 
must be positive in the same z region as shown by the dotted curve. The 
distribution Ag(0.05/z) in the region 0.05 < z < 0.1 corresponds to Ag(x) 
in the region 0.5 < x < 1. Therefore, the gluon distribution should be 
positive at large x for fitting to the experimental data. 

3. Summary 

For determination of A g(x), we performed the global analyses with present 
DIS and ir° production data. Although the uncertainty of the fist moment 
is significantly reduced by adding the tt° data, the sign problem of Ag(x) 
appears. However, the DIS data of COMPASS and HERMES experiments 
provide a constraint on the large- a; behavior of Ag(x) through the NLO 
gluon correction term in g\ (x) . 
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